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Summary.— The deep -inelastic deuteron structure function (SF) F^xd.Q 2 ) in 
the covariant approach in light-cone variables is considered. Thexp and Q 2 -dependences 
of SF are calculated. The QCD analysis of generated data both for non- cumulative 
xd < 1 and cumulative x^ > 1 ranges was performed. It was shown that Q 2 -evolution 
of SF is valid for ranges 0.275 < xd < 0.85 and 1.1 < xd < 1.4 for the same value 
of QCD scale parameter A. It was found the xd- dependence of SF for the ranges is 
essentially different. 

PACS: 12.38. Qk; 12.38.Bx; 12.90. +b. 



1.— Introduction 

The progress of perturbative Quantum Chromodynamics (QCD) in the description 
of the high energy physics of strong interactions is considerable jl]. Deep- inelastic 
scattering (DIS) of leptons provides the data on nucleon structure functions (SF) which 
could be used for the most precise comparison with QCD. The accuracy of SF data is 
sufficient to be sensitive to nuclear effects in DIS, e.g. the EMC effect [0. A direct 
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evidence of the nucleus structure in DIS has been obtained by measuring the carbon 
SF [H] in the kinematical region forbidden for scattering on a free nucleon (the Bjorken 
kinematical variable x = Q 2 /2(pq) is more than 1 in the range). The region is known 
as a cumulative region |4], |5| and it is intensively investigated both experimentally 
and theoretically in proton- nucleus and nucleus-nucleus collisions [(J . An important 
matter of the relativistic nuclear physics is to adduce arguments for applying QCD to 
analyse SF at xp > 1. The deuteron as a simplest nuclear system is an excellent object 
to investigate the problem of relativistic theory of a deuteron. 

Dependence of the deuteron structure functions _F 2 D onic in the cumulative range 
has been studied by many authors. Different models (the few-nucleon correlation model 
0, the ffucton model || ) are used to describe a high momentum component of the 
deuteron wave function. The microscopic picture for cumulative phenomena on the 
deuteron near the threshold based on perturbative QCD has been developed in ][|. 

The dependence of the nuclear structure function ratio Rpl D = F 2 C ' / F® on xd 



has been considered in [10]. In the analysis the BCDMS experimental data on F : 
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H have been used. The deuteron SF has been calculated in the light-cone covariant 
approach [jlT| and it was found that the ratio R c / D {xo) reveals an exponent growth 
~ exp(a-xr>) for xp > 1 with the slope parameter a ~ 6.6 and differs significantly from 
the ratio behaviour for xd < 1. The latter is in good agreement with experimental 
data F3, [TJ]]. It is considered that the ratio R c / D for xjj < 1 describes quark 



distributions in a nucleon of a nucleus. The ratio R C I D for xd > 1 describes quark 
distributions in a nucleus. In JTH] similar exponent dependence of the inclusive pion 
backward cross section ratio Rrl D ~ p\(X)/pJ ) (X) on the cumulative number X 
|16| for the p + A — > 7r(180°) + ... process was found. It was concluded that this 



dependence demonstrates the change of the regime from hadron to quark degrees of 
freedom. We can consider that the exponent growth of the similar ratio R A I D for 
different processes (deep-inelastic lepton scattering, deep-inelastic nucler reactions) in 
the cumulative range is the general feature of quark interactions in the superdensity 
matter. 

An important matter in the description of DIS is whether or not the Q 2 -evolution 
of the structure function is possible in the framework of QCD in the cumulative range. 
If it is possible, the QCD analysis could be used to study the transition regime and to 
determine quark distribution functions in nuclei. 

In the present paper, the QCD analysis of data on the deep-inelastic deuteron struc- 
ture function F®(xd,Q 2 ) for xd < 1 and xd > 1 is performed. The data on the xd 



and Q -dependence of SF are simulated by the relativistic deuteron model |TT[. The 



simulated data for F®(xd,Q 2 ) for the range xp < 1 are in good agreement with ex- 
perimental data WM- The experimental data for the deuteron SF in the range xd > 1 
for high Q 2 are absent. These data would be very interesting both to verify deuteron 
models and methods of describing relativistic nuclear systems and to study the high 
momentum component of the relativistic deuteron wave function (RDWF). The QCD 
analysis of data is based on the expansion method of SF over Jacobi polynomials ||17|| - 
JTJJ. The Q 2 -evolution of the moments of the structure function F® is found as the 
solution of the corresponding renormalization equation. Effectivity of the method to 
control higher perturbative QCD corrections and to investigate sensitivity of the QCD 
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scale parameter A were shown in [ljj. It was found that the data on F^(xd, Q 2 ) obey 



the (^-evolution of SF in the framework of QCD both for 0.275 < xd < 0.85 and 
1.1 < xd < 1-4 ranges in accordance with the criterium \ 2 very well. The value of the 
parameter A is determined to be the same for both the regions. We would like to note 
that the parameters of structure function parametrisations differ significantly for these 
ranges and this gives us some evidence for two regimes of the ^-dependence of the 
structure function (xd,Q 2 ). 

2.— Deuteron Structure Function F^ixD-, Q 2 ) 

The deuteron structure functions Wi^i^iQ 2 ) are related to the imaginary part of 
the forward scattering amplitude of a virtual photon on a deuteron by the standard 
formula 



(1) W£ = -{g»„ - q»qv/q 2 ) ■ W? + (p„ - q,(pq) /q 2 )( P u - qu{pq)/q 2 ) • W 2 D /M 2 . 

Here q, p are momenta of a photon and a deuteron; M is the deuteron mass. 

In the relativistic impulse approximation (RIA) the forward scattering amplitude of 
the virtual photon 7* on the deuteron A® is defined via a similar scattering amplitude 
on the nucleon A^ u as follows 

(2) A°{q,p) = J ^SpiA^kJ-T^h)}. 

In the expression (2) T(si, ki) is the amplitude of the forward N — D scattering and 
the conventional notations Q 2 = —q 2 > 0, v = (pq), Si = {p — k\) 2 are used. The 
integration is carried out over the active nucleon momentum k\. The calculation of the 
imaginary part of the amplitude A® gives us the possibility to put the nucleon spectator 
with the momentum k = p — k\ on mass shell. Therefore the tensor W^ u is expressed 
via the DNN vertex with one nucleon on mass shell. The vertex is described by the 
function T a (ki) and depends on one variable k±. The vector index a characterizes the 
deuteron spin. With the relation between the RDWF and the vertex function T a (ki): 
i^aiki) = T a (ki) ■ (m + k\)~ l , the expression for the tensor can be written as 



= J ^5(m 2 -k 2 )Q{k )Q(p + -k + )Sp{w N ^ 

Here Pap(p) = —(g a /3—PaPa/M 2 )/3 is the deuteron polarization density matrix and the 
- function and light-cone variables (k± = k Q ± fc 3 , k±) are used. The vertex function 
T a (ki) is defined via 4 scalar functions a,i{k\) (i — 1 — 4) and has the form 

(4) r a (k!) = k la [ ai {kl) + a 2 (^)(m + fcj)] + 7a[a 3 (^ 2 ) + a 4 (^)(™ + fci)]. 



The relativization procedure of deuteron wave function ip a has been proposed in [TT, pO 
The scalar functions a>i(k\) have been constructed in the form of a sum of pole terms. 
Some pole positions and residues have been found from the comparison of our RDWF 
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in the nonrelativistic limit with the known nonrelativistic deuteron wave function. 



For the latter the Paris wave function [21] was taken. Other parameters were fixed 



from the description of the static characteristics of the deuteron (an electric charge - 
G e (0) = 1(e), magnetic - G m (0) = /ic(e/2Mo) and quadrupole - Gq(0) = Qu^/Mp) 
moments) in the relativistic impulse approximation. 

The calculation of (3) in the light-cone variables gives the final expression for the 
deuteron SF F 2 D = vWf 

(5) Fj>(x D , Q 2 ) = J 1 dx d 2 k ± p(x, k ± ) ■ F?(z/x, Q 2 ), 

where xp = 2z, < z < 1. The nucleon SF = (F$ + F£)/2 is defined by the 
proton and neutron ones. The positive function p(x, kj_) describes the probability that 
the active nucleon carries away the fraction of the deuteron momentum x = ki + /p + 
and the transverse momentum k± in the infinite momentum frame. It is expressed via 
the vertex function r a (A;i) as follows 

(6) p(x, k ± ) oc Spffifa) ■ (m + k) ■ ^ p {k x ) ■ q ■ p a p{p)}- 

Note that in the approach used the distribution function p(x, k±) includes not only 
the usual S- and Z)-wave components of the deuteron but the P-component too. The 
latter describes the contribution of the iViV-pair production. 

In the RIA the deuteron SF F 2 D is defined by (5) as a sum of the proton and neutron 
SF integrated on the x and k±. The NMC data @ on the ratio Rp /p = F<?/F% and 
F% and the relativistic deuteron model have been used to extract the neutron SF F£. 
For the latter the parametrisation 

(7) F 2 n (x, Q 2 ) = (1 - 0.75x)(l - 0.15^(1 - a;)) • F 2 p (x, Q 2 ) 

has been obtained in [Q. The results for the absolute value of _F 2 D calculated within 
the parametrisation of the neutron SF are in good agreement with experimental data 
both in low |TjJ and high [12, |T^| x/j-ranges. 



3.— Proton Structure Function FP(x : Q 2 ) 



In our analysis we shall use the parametrisation of F^x^Q 2 ) given in [IB]]. The 
parametrization describes the NMC, SLAC and BCDMS data very well. The verifi- 
cation of this fit in the region 0.006 < x < 0.6 gives A = 200 MeV. To parametrize 
F2(x,Q 2 ) at 0.55 < x < 1, we have made the QCD fit of parametrisation [13| con- 
sidering it as an "experimental points" at 0.275 < x < 0.55 with the leading order 
nonsinglet evolution of the moments of F%{x, Q 2 ) with A = 200 MeV . The SF at the 
fixed point Ql = 10 (GeV/c) 2 was parametrized as 

(8) Ff (x, Ql) = Ax B (l - xf (1 + 7 z). 

The parameters A, B, C and 7 in Eq. ([8]) are free parameters and are determined 
by the fit of the data. Then, on the basis of expression (§) the values of SF for 
0.55 < x < 1 were calculated. To achieve agreement between QCD - evolution results 
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and parametrization the former should be multiplied by the factor R(x,Q 2 ) for 
0.55 < x < 1: 



(9) R(x, Q 2 ) = (1. + 17.611(x - 0.55)- (3 - +0 - 661 ^ (Q /Q o))) (x - 0.55). 

The modified parametrization of the proton SF in a wide range of the Bjorken 
variable x is used as an "experimental" input for the deuteron model described above 
to simulate the data on the deuteron SF. The data for the deuteron SF are simulated 
for x D = 0.l- 1.8 and Q 2 = 17 - 230 {GeV/cf. 

Figure 1 shows the x^-dependence of F 2 d (xd, Q 2 ) at Q 2 = 61.5 (GeV/c) 2 . One 
can see that the SF falls down drastically for xd > 1- The open points are BCDMS 
data for SF F 2 for the /i + 12 C — > \i + ... process [|J. The model-simulated data for 
F 2 D (xo, Q 2 ) are in a qualitative agreement with BCDMS data. We would like to note 
that the more detail calculations |TD[ of the dependence of the nuclear structure func- 



tion ratio R c/D = F^/F^ onic with RDWF [|ll| have shown that the ratio is similar 
to the experimental data for 0.01 < xd < 1 P3[ and exponentially growthes for xd > 1- 



4.— Method and Results of Analysis 

Now we can apply the nonsinglet QCD fit to the data on the deuteron SF simulated 
in the previous section. We have reduced the deuteron data to a standard interval of 
the scaling variable < x < 1 by putting x = xd/2 and have applied the method of 
the QCD analysis based on the expansion of SF over the Jacobi polynomials |p"Tf-||19|. 
The SF is presented as follows 

*V max 1 

(10) ff— (x } Q 2 )=x a (l-xf £ Q^x)Y,c < f\a J (3)Mj + 2(Q 2 ),n = 2,3 } ... 

n=0 j=0 



II (<?- - asiQ2o) 



M n (Ql) 



_oc s (Q 2 

where d n = — / y% s /2f3 , (3q = 11 — §-/V/, Nf is the number of flavours, 7^5 is the 
anomalous dimension in the one-loop perturbative QCD approximation: 
7 $j = 8/3 [4S-i(7i) - 3 - 2/(n(n + 1))], where S^n) = £?=i 
The moments M n (Q 2 ) of SF are defined by 

1 



(12) M n (Q 2 ) = I dx x N ~ 2 F 2 D (x,Q 2 ). 

Jo 

The parametrization form of SF similar to (8) is used at a fixed point 



(13) F 2 D ( Xj Ql) = A(x D /2f(l - x D /2f (1 + jx D /2). 

The constants A, B, C, 7 in ( |I3|) and the QCD scale parameter A are considered as 
free parameters to be determined at Qq. 

The results of the fit are presented in Table 1. We have found that values of A 
are approximately the same in two different regions of xd and not significantly differ 
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from the input proton parameter A = 200 GeV ||13|| . The xd- dependence of SF for 
two ranges xd = 0.225 — 0.85 and xd = 1.1 — 1-4 , is determined separately. It is 
established that the shapes of SF are essential different (see Table 1). The obtained 
value of Xd.f. corresponds to 10% for errors of "experimental" deuteron data. 



Table 1. The results of the LO nonsinglet QCD fit of the deuteron SF. Nf = 4, 
17 < Q 2 < 230 (GeV/c) 2 , N max = 12. The errors of F 2 D are put to be 10%. 





0.275 <x D < 0.85 


1.1 <x D < 1.4 


A [MeV] 


268 ± 67 


240 ± 44 


A 


8.39 ± 1.50 


(6.26 ±0.84)10~ 5 


B 


0.898 ± 0.065 


-0.115 ± 0.015 


C 


8.659 ± 0.456 


11.07 ± 0.07 


7 


-2.034 ± 0.055 


3376. ± 457. 


Xd.f. 


12.0/121 


34.4/68 



Figure 2 shows the dependence of F®(xd, Q 2 ) on Q 2 for xd = 1.1 — 1.8. The results 
of the QCD fit for xu = 1.1 — 1.4 are drawn by solid lines, and the results of direct 
calculation by formula (5) are shown by open points. Good agreement between the 
model-simulated data on F^(xd, Q 2 ) and results of the QCD fit is observed. 

We would like to note that a more general task is a simultaneous QCD fit of _F 2 D in 
the whole range xd = 0. — 2. but it requires a rather complicated a;-parametrization of 
SF instead of (8). 

5.— Conclusion 

The QCD analysis of data on the deep-inelastic deuteron structure function F®(xd, Q 2 ) 
was performed. The data were simulated by the relativistic deuteron model in the co- 
variant approach in light-cone variables. It was found that the data on F^(x£i,Q 2 ) 
obey the Q 2 -evolution of SF in the framework of QCD both for 0.275 < xo < 0.85 
and 1.1 < xd < 1-4 in accordance with the criterium \ 2 very well. The value of the 
parameter A was determined to be the same for both the regions. The parameters of 
structure function parametrizations were determined and it was found that they differ 
significantly for these ranges. These results give some evidence for two regimes of the 
a; ^-dependence of the structure function F®(x[),Q 2 ) . 

We would like to emphasize that the QCD Q 2 - evolution of SF in the cumulative 
range (xd > 1) could be a crucial test for verification of both the nuclear models and 
the parton model, and thus measurements of SF in DIS for xd > 1 should be performed 
with a high accuracy at CERN, DESY and FERMILAB. 
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Figure captions. 



Fig 1. Deep-inelastic deuteron structure function F®(xd,Q 2 )- Theoretical results 
have been obtained with the proton parametrization F$ taken from [13]. Experimental 
data for F 2 C are taken from [3]: o — 61.5 (GeV/c) 2 . 

Fig 2. Deep-inelastic deuteron structure function F®(xd,Q 2 )- Open points are 
generated "experimental" data; solid lines are results of the QCD fit of the data. 
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